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In  this  paper  we  will  give  a  simple  approach  to  determining  conditions 
for  stability  of  linear  feedback  systems  subject  to  additive  and  multipli¬ 
cative  perturbations  in  the  operators  describing  these  systems.  The  approach 
is  based  on  techniques  used  in  functional  analysis,  and  provides  an  alternative 
development  and  generalization  of  some  conditions  for  the  time  -  invariant 
case  that  have  appeared  in  the  literature  very  recently.  As  an  example  of 
the  application  of  the  conditions,  we  consider  the  determination  of  finite 
regions  of  stability  for  singularly  perturbed  systems.. 


1.  Introduction 


An  important  theme  in  system  theory  is  the  preservation  of  various 
system  theoretic  properties  in  the  face  of  variations  in  the  system  model. 
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It  is  possible  tc  distinguish  twc  variations  on  this  theme.  In  the  first, 
attention  is  restricted  to  infinitesimal  changes  in  the  parameters  of  the 
nominal  system  model.  Thus  cne  begins  by  assuming  that  the  nominal  system 
nas  a  certain  property,  and  then  asks  if  there  exists  an  open  set  about  the 
nominal  system  parameters  such  that  all  the  s'/stens  with  parameters  in  this 
set  have  the  desired  property.  V-  will  refer  to  investigations  of  tnis  first 
type  as  sensitivity  theory.  A  second  approach  requires  the  explicit  delmeati 
of  finite  regions  of  models  about  the  nominal  model  fcr  which  the  given 
property  is  preserved.  We  will  refer  tc  investigations  of  this  second  type 
as  robustness  theory . 

Within  the  context  of  sensitivity  or  robustness  theory,  there 
are  several  properties  that  have  been  investigated,  For  example,  it  is 
well  known  that  the  controllability  property  is  insensitive  to  small  para¬ 
meter  variations  II, p,  43] .  As  another  example,  it  is  well  known  that 
type-1  servomechanisms  have  zero  steady  state  step  tracking  error  despite 
large  (but  not  destabilizing)  variations  in  their  transfer  function  matrices 
[2]  . 

In  this  paper  we  will  focus  on  the  robustness  of  the  stability 
property  of  linear  multivariable  feedback  systems.  This  subject  is  of  special 
interest,  since  stability  is  the  most  basic  system  theoretic  issue  and  since 
practical  feedback  systems  must  remain  stable  in  the  face  of  large  para¬ 
meter  variations. 

The  importance  of  obtaining  robustly  stable  feedback  control 
systems  has  long  been  recognized  by  designers  [3] .  Indeed,  a  principal 
reason  for  using  feedback  rather  than  open-loop  control  is  the  presence  of 
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model  uncertainties.  Any  model  is  at  best  an  approximation  of  reality,  and 
the  relatively  low  order,  linear  time-invariant  models  most  often  used  for 
controller  synthesis  are  bound  to  be  rather  crude  approximations. 

In  classical  frequency  domain  techniques  for  single-input,  single-out¬ 
put  (SISO)  control  system  design,  the  robustness  issue  is  naturally  handled  (3]. 
These  techniques  employ  various  graphical  means  (e.g.,  Bode,  Nyquist,  inverse- 

Nyquist,  Nichols  plots)  of  displaying  the  system  model  in  terms  of  its  frequency 
response.  From  these  plots,  one  can  determine  by  inspection  the  minimum  charge 
in  the  model  frequency  response  that  leads  to  instability.  These  changes  are 
often  quantified  by  the  gain  and  phase  margins  of  the  feedback  system;  sometimes 
the  design  is  required  to  have  certain  minimal  margins  in  order  to  be  acceptable 
[4,  p.  43]. 

In  modern  time  domain  techniques  (such  as  the  pole  placement  or 
linear-quadratic-Gaussian  approaches)  for  multiple-input,  multiple-output 
(MIMO)  systems,  the  robust  ess  issue  is  not  directly  dealt  with.  Instead,  it 
is  necessary  to  transform  the  resulting  design  to  the  frequency  domain  to  examine 
its  robustness  properties.  For  SISO  systems,  this  is  accomplished  as  for 
classical  designs,  but  the  situation  is  less  clear  in  the  MIMO  case,  where  it 
is  necessary  to  consider  simultaneous  variations  in  the  frequency  responses 
of  all  the  loops. 

Very  recently,  there  has  been  some  important  work  addressing  the 
multivariable  robustness  issue.  In  his  thesis  [5,6]  Safbnov  gives  a  powerful 
approach,  based  on  a  multivariable  sector  stability  theorem,  that  can  characterize 
robustness  for  very  general  nonlinear  MIMO  feedback  systems.  In  a  recent  paper 
[7] ,  Doyle  develops  a  robustness  characterization  for  the  linear  time  invariant 
MIMO  case  (It  can  be  shown  that  Doyle's  result  can  also  be  obtained  by  Safonov's 
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approach  [8].)  Doyle's  characterization  involves  computing  the  minimum  singular 
value  of  a  certain  transfer  function  matrix,  and  this  computation  essentially 
determines  the  minimum  simultaneous  variation  of  the  system  frequency  responses 
that  leads  to  instability .  Since  there  is  sophisticated  and  widely  accessible 
software  to  compute  singular  values  [9],  this  characterization  is  of  great 
practical  value. 

The  present  paper  is  prompted  by  two  observations.  First,  the  use 
of  singular  values  to  characterize  robustness  is  suggestive  of  connections  with 
numerical  analysis,  but  these  connections  are  not  clear  from  Doyle's  approach 
utilizing  the  multivariable  Nyquist  theorem.  Second,  a  specific  instance  of  the 
robustness  question  arises  when  a  system  is  approximated  by  making  a  singular 
perturbation  to  reduce  its  order,  A  related  motivation,  although  only  briefly 
discussed  in  this  paper,  arises  from  the  desire  to  use  multi-model  techniques  in 
the  design  of  decentralized  controllers  for  large  scale  systems  [10,11]. 

The  structure  of  this  paper  is  as  follows.  In  Section  2  we  consider 
the  robust  stability  of  MIMO  linear  feedback  systems  using  a  generalized  numerical- 

analytic  approach.  When  specialized  to  the  time  invariant  case,  with 
rational  transfer  function  matrix  perturbations,  Doyle's  characterization 

results.  In  Section  3  we  will  apply  the  results  of  Section  2  to  a  specific 
robustness  question  arising  in  singular  perturbation  theory.  Section  4  contains 
the  summary  and  conclusions. 

Notation 

We  will  use  the  standard  notation  of  input-output  itability  theory: 
see  (12,  pp.  13-14],  or  [13,  pp.  38-39], 

X  «  some  Banach  space  of  functions  x:T  -*•  X 
T  »  subset  of  the  real  numbers 
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X  =  finite  dimensional  vector  space 
X^  =  {x  :  ?_  x  e  X  for  an  tst| 

f  x  ( t) 


(PTx)(t)  = 


{ 

10 


t  <_  T 
t  >  T 


=  space  of  m- vector  functions  on  T 
with  integrable  Euclidean  norm 


I  :  X^  ■»  Xe  =  identity  operator 


G  : 


is  causal  if  P_GP_  =  P_G 
for  al)  T 


T  T 


A*  =  conjugate  transpose  of  a  complex 
matrix  A. 


2.  Robust  Stability  of  Linear  Systems 

We  consider  the  feedback  system  depicted  in  Figure  1.  Here 

the  causal  linear  operator  G  :  S^e  represents  the  plant 

plus  any  compensation  that  is  used.  The  basic  feedback  equation 
is 

(I  +  G)e  -  u  (2.1) 

and  the  basic  stability  question  is  whether  (I  +  G)  1  "*'i?2e 

exists,  is  causal,  and  is  a  bounded  linear  operator  when  restricted  to  the 
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subspace  2^  of  2^  .  We  will  assisie  that  the  nominal  system  is  stable  in 

n 

this  sense  throughout  this  section,  l.a.  for  u  e  there  exists  a  uniaue 

2e 

causally  related  e  £  2^^  satisfying  (2.1),  and  that  u  implies  that  the 

corresponding  e  and  consequently  y  =  u-e  z2^.  We  are  interested  in 

whether  the  closed  loop  system  retains  these  properties  when  subject  to 
additive  (Figure  2)  or  multiplicative  (Figure  3)  perturbations  representing 
uncertainty  in  the  dynamical  behavior  of  the  system. 

The  following  theorem  provides  the  basis  for  our  analysis. 


THEOREM  1 

Let  A:  -*■  be  a  causal  linear  operator,  and  suppose  A  ^ 

exists,  is  causal,  and  is  bounded  when  restricted  to  X  .  Then,  if 

AA  :  X  -►  X  is  a  causal  linear  operator  that  is  bounded  when  restricted 
e  e 

to  X  ,  and  if 

| | A_1AA |  |  ^  <  1,  (2.2) 

it  follows  that  (A  +  AA)  ^  :  X  -►  X  exists,  is  causal  and  is  bounded  when 

e  e 

restricted  to  X . 

Proof 

The  operator  A  ^AA  is  well  defined, causal,  and  bounded  on  X 
by  assumption.  Since  | | A  ^AA j | ^  <  l,  the  contraction  mapping  theorem 

implies  that  the  sequence  x^,  k  =0,1...,  defined  by 

Vl  *  -A^AAx*  +  b;  xo=0  (2.3) 


converges  to  a  unique  solution  xeX  of 
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J  +  A_1lA)x 


:  2 . 4 ) 


for  any  b  s  X.  Thus i I  *  A  ^£Aj  ^  :  X  -  X  exists,  and  is  therefore  bounded 
since  I  +  A  ’''AA  is.  Causality  follows  since  each  iterate  depends  causa 
on  b  and  consequently  the  limit  x  of  x^  depends  causally  on  b.  Since 
(I  +  A  \lA)  ^  is  causal,  it  can  be  uniquely  extended  to  X„  by  requiring 


P  (I  +  A~1AA)“1x  =  ?  (I  4  A_1AA)_1F  x 

T  T  t 


(2.5) 


for  x  £  Xg.  Finally,  defini 


ng 


(A  +AA)"1  =  (I  +  A"11A)*'1A"1, 


(2.6) 


gives  the  required  inverse  of  A  +  AA.  Q.E.D. 


Remarks 


1.  Since 

I  I  A_1AA|  |  x  1  i  i A_1 !  ix  i!Mi!x  (2.7) 

a  sufficient  condition  for  (2.1)  is 

i  j  A"* 1 1  |x  [  I AA I  I  x  <  1.  (2.8) 

2.  A  basic  result  in  numerical  analysis  is  that  if  an  nxm  matrix  A 
is  invertible,  then  A  +  AA  is  invertible  for  all  AA  satisfying 

a(AA)  <  a  (A)  (2.9) 

where  a  (AA)  =  ||AA||2,  a  (A)  =  ( i  i  a”1;  j  )_1 

are  respectively  the  smallest  and  largest  singular  values  of  the  matrix  A. 
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Theorem  1  is  thus  a  generalization  of  this  classical  finite  dimensional  result  wher 
of  course,  boundedness  and  causality  are  not  at  issue.  In  the  finite  dimensional 
case  there  exists  AA  such  that  z  (AA)  =  c_tA)  and  A  +  iA  is  singular;  this  is 
easily  proved  by  the  singular  value  decomposition^. 

3.  The  contraction  mapping  argument  in  Theorem  1  is  a  standard  technique  of 

applied  mathematics ;  the  fact  that  X  is  not  a  Banach  space  complicates  the  argumen 

e 

The  causality  argument  has  been  used  by  Willems  [12,  p.  93]  in  a  slightly 
different  context.  The  linearity  of  the  perturbation  operator  is  not  essential. 

4.  Theorem  1  can  be  used  to  give  SC  robust  stability  results,  but  we  will 

P 

confine  ouselves  to  the  case  p  =  2  in  the  sequel. 

5.  Theorem  1  can  be  used  to  obtain  robust  stability  results  for  both 
continuous  and  discrete  time,  but  in  the  sequel  we  will  confine  our  attention 
to  the  case  T  ■  [0 ]  . 

The  robust  stability  questions  posed  at  the  beginning 
of  this  section  are  now  answered  in  terms  of  Theorem  2. 

THEOREM  2 

Assume  that  the  basic  feedback  system  of  Figure  1  is  stable.  Then 
(i)  the  system  remains  stable  for  additive  pertubations  G  (Figure  2)  provided 

|  |  (I  +  G)_1AG|  '  m  <  1  (2.10) 

"2 

and  (ii)  the  system  remains  stable  for  multiplicative  perturbations  AG  (Figure  3) 

^The  singular  value  decomposition  of  an  nxn  nonsingular  complex  matrix  A  is 

A=UE V* ,  where  U  and  V  are  unitary  nxn  matrices,  E  *  diag  (a  ,  . ,0  )  and 

the  singular  values  a.  are  the  non-negative  square  roots  of  the  eigenva?ues  of 
A*A.  See  [14]  for  references,  a  more  general  definition,  and  an  excellent 
discussion  of  the  fundamental  role  of  the  singular  value  decomposition  in  linear 
systems  theory. 
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provided 

i  !  (I  -  (I  +  G)“1]1G|  \(j*  <  l  (2.11) 

Proof 

For  case  (i) ,  we  apply  Theorem  i  to  the  equation 

(7  +  G  +  AG) e  =  u  (2.12) 

while  for  case  (ii)  we  consider 

(7+G(7+AG)Je  =  ( I+G+ ( 1+0) AG-AG) e  =  u.  Q.E.D.  (2.13) 

The  practical  importance  of  Theorem  2  stems  from  the  fact  that  the 
2^  norm  of  a  linear  convolution  operator  can  be  computed  from  its  transfer 
function  matrix.  This  fact,  which  is  a  consequence  of  Parseval’s  Theorem, 
is  well  known  in  the  input-output  stability  theoretic  literature;  see, 
e.g. ,  [13,  p.  26]  . 


.eirjna  1 


Let  the  operator  G;  ■*  2^  for  T  =  [0,»]  be  define!  by 


(Gx)  (t 


G ( t-T  )  X (T  )  dT 


(2.14) 


where  the  elements  of  the  impulse  response  matrix  G(t)  are  assumed 
absoiutely  integrable  on  T.  Then 


where 


a 


a 

max 


OilGO)) 


(2.15) 


max 


max  max 

<jj>0  ldi^.m 


(2.16) 
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and  where  0^£G(j'i))  denotes  the  ith  singular  value  of  the  transfer  function 
matrix  corresponding  to  G  . 

Combining  Theorem  2  and  Lemma  1,  we  obtain  the  following  result. 
THEOREM  3 

Assume  that  the  nominal  system  in  Figure  1  is  time  invariant,  stable, 
and  that  the  operators  (I+G)  ^ ,AG  can  be  represented  as  convolution 
operators  with  impulse  response  matrices  with  absolutely  integrable  elements. 
Then  (i)  the  system  remains  stable  for  additive  perturbations  “G  satisfying 

o’  (A  G  (  jco )  )  <  £(I+G(jc»)),  b  >  0  (2.17) 

and  (ii)  the  system  remains  stable  for  multiplicative  perturbations 
satisfying 

o(AG(  jco)  )<  a(l+G_1(jaj) )  ,  oj  >  o.  (2.18) 

Here  G(jaj)  and  AG(jw)  are  the  transfer  functions  of  G  and  AG,  and 
0(A)  and  0(A)  denote  the  maximum  and  minimum  singular  values  of  A. 

Proof 

(i)  From  Lemma  1  and  (2.17)  we  nave 

| |AG| I  <  II  (i+G)  L| ,  1  (2.19) 

so  that 

|  |  (nG)"1]  |  |  |aG|  |  <  1.  (2.20) 

(ii)  Note  that 

a  (I+G-1  ( ja>) )  =  |  a  f(l+c-_1(  jco)  )_1]  }  1 

-  |  ff[(i+G(jcj))  _i  G(  joj)  ]  I’1  =  a[l-(X+G(ja)))"1]  (2.21) 

Therefore  Lemma  1  and  (2.18)  imply  that 

I  I  I  -  (I+G)_1|  j  I |AG| I  <  1.  (2.22) 


Q.E.D. 


Remarks 


1.  Notice  the  analogy  between  Theorem  3  concerning  robust  stabilitv  ci 
linear  systems  and  the  classical  result  quote!  nreviousiy  concerning  robust 
inversion  of  matrices  (or  bounded  operators!. 

2.  Theorem  3  for  the  case  of  rational  transfer  function  matrices  is  me 
result  of  Doyle  alluded  to  previously.  Doyle's  proof  is  completely  different, 
aowever,  depending  on  cne  multivariable  Nyquist  Theorem,  so  that  the  connect: 
wj  tli  the  inversion  issue  are  only  implicit. 

3.  The  quantity  £_(  I+G(  jb  )  )  or  3  (1+3  ^  1 3 -h  )  is  easily  computed  and  plctte 
as  a  function  of  gj.  Doyle  has  made  great  use  of  this  technique  in  the  analys 
of  multivariable  feedback  systems.  Such  a  nlot  plays  much  the  same 

role  for  determining  MIMO  robustness  properties  as  the  more  classical  Bode, 
etc.  plots  in  SISO  design. 

4.  The  quantity  a  1  ! I+G ( jco)  )  is  the  generalization  of  the  classical  Bode 
SISO  sensitivity  function  of  changes  in  the  closed-loop  transfer  function  with 
respect  to  changes  in  the  open- loop  transfer  function  in  the  following  sense. 
Let  y^  denote  the  output  of  the  system  of  Figure  1  for  a  given  input  and 

y ^  the  corresponding  output  of  the  system  of  Figure  2  for  the  same  input. 

Then  one  can  show 

y1(jcu)  -  y2<  jcu)  *  (I+G(  jcu)  (G(jo>)  +  AG(  joj) )  ^(jcu) 

(2.23) 


so  that 

Jy^jco)-  y2(jw)||<  ;;  «.G(  jw)  (l+G(jo»)  +  AGCjcu)  )_1i  u 

i  1  y 2  <  jou )  ’,  I  •  (2,24) 

Consequently,  the  percentage  change  in  the  closed-loop  transfer  function 
matrix  i3  attenuated  from  the  percentage  change  in  the  open  loop  transfer 


function  matrix  by  the  factor  o  *(:  -  G(j._t) 

5.  On  the  other  hand,  the  perturbations  defining  the  SISO  gain  and 
phase  margins  are  multiplicative  rather  than  additive,  so  that  7(1  +  i  A, 
is  more  appropriate  as  a  measure  of  the  tolerance  of  the  feedback,  system 
to  model  uncertainty. 

6.  It  is  in  general  impossible  to  express  a  (I  +  G  ( jiu )  )  in  terms  of 
a (  I  +  G(  jaj)L 

7.  Consider  the  feedback  system  in  Figure  1  with 

G(s)  =  -G(sl  -  A)_1B  (2.25) 


where 

G  =  B'K  (2.26) 

0  =  A'K  +  KA  +  C'C  -  KBB'K  (2.27) 

(We  assume  [A,B]  controllable  and  [A,C]  observable  so  that  a  unique  positive 
definite  solution  of  (2.27)  exists.)  The  well  known  equality  [15] 

[I  +  G(-sI  -  A)_1B]'  [I  +  G(sl  -  A)_1B)  *  I  +  IC (-Sl-A)_1B]'  [C(sl  -  A)_1B], 

(2.28) 


which  follows  from  (2.27)  after  a  little  manipulation,  shows  that  the  system 
is  robust  to  additive  pertubations .  Safonov  and  Athans  [6]  have  shown  that 
the  system  of  Figure  3  with  G(s)  defined  by(2.25)is  stable  for  AG(jcj)  satisfying 
0(AG(jco))  £  1/2  for  all  co  _>0; 

this  can  also  be  inferred  directly  using  the  inequality  [19J 


a(I  + 


G  1  ( ja>) ) 


0(1  +  G(  joj)  ) 

1  +0(1  +  G(  jco)  ) 


(2.29) 


together  with  (2.28)  and  Theorem  3 


Application  to  Singular  Per turbaticr.  Theory 


3  . 

In  the  previous  section  we  have  discussed  the  robustness  of  the  stability 
property  of  a  linear  dynamic  system,  to  model  variations.  In  this  se'tnr.  -*v 
will  consider  a  particular  form  of  mcdel  variation  due  to  a  singular  perturrati: 
We  consider  systems  of  the  form 


xl ;  t)  =  +  Ap-.x2  c) 

£X2(t)  =  A21X1U)  +  A22X2lt) 


(3.11 


where  e?0  is  a  small  parameter,  and  it  is  assume  that  the  matrix  A  ^  is 
stable  (has  eigenvalues  with  negative  real  parts) .  We  define  the  so  called 
degenerate  system 


(AU  '  A12A22  A21)Xld(t) 


(3.2) 


associated  with  (3.1).  This  system  is  a  reduced  order  system  that  neglects 
certain  high-frequency  or  parasitic  effects  incorporated  in  the  model  (3.1). 

It  has  been  shown  that  the  stability  of  (3.2)  (in  the  sense  that  the  eigen¬ 
values  of  the  system  matrix  have  negative  real  parts)  is  insensitive  to  these 

effects  in  the  sense  that  there  exists  £  >0  such  that  (3.1)  is  stable  for 

o 

all  0<£<£  if  (3.2)  is  [17].  We  propose  to  examine  the  robustness  of  the 
o 

stability  (in  the  input-output  sense  of  Section  2)  of  (3.2)  to  the  parasitic 
effects  present  in  (3.1)  . 

We  begin  by  Laplace  -  transforming  equations  (3.1)  (assuming  zero  initial 
condition) . 


1 


See  [16]  for  an  excellent  survey  of  results  in  singular  perturbation  theory 
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To  apply  the  input-output  stability  results  of  the  previous  section  it  is 
necessary  to  apply  a  test  input  to  the  system.  This  is  most  conveniently 
done  as  illustrated  in  Figure  4,  although  other  locations  are  possible^. 

Figure  4  closely  resembles  Figure  3  with 


G(s) 


A21(S1 


(3.5) 


AG(s,£)  =  — £s  ( £ S I  -  A  ^ 


(3.6) 


except  that  the  perturbation  is  post-multiplicative  rather  than  pre- 
multiplicative.  However,  assuming  G(s)  has  full  rank  as  a  rational  matrix, 
it  is  easily  verified  tnat  the  analysis  of  the  preceding  section  is  essentially 
unaffected.  Thus  we  have  the  following  result. 


To  insure  the  equivalence  of  the  input-output  stability  analysis  with  the 
condition  that  the  system  matrix  of  (3.1)  has  eigenvalues  with  negative  real 
parts,  it  is  necessary  to  have  the  conditions 
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controllable , 
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Figure  4. 


Singular  Perturbation  in  the  Frequency  Domain. 
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THEOREM  4 

Assume  that  the  system  of  Figure  4  is  stable  for  e  =  0.  Then  it  remains 
stable  for  all  e  >  0  satisfying  the  inequality 


a  (AG ( joj,e) )  <_  a(l  +  G-1(  jto ) ) 

for  all  co>0. 

The  use  of  Theorem  4  is  illustrated  by  the  following  examples. 


(3.7) 


Example  1 
x^  (t) 

ex  (t) 

G(s)  = 
AG(s,e) 


x  (t)  +  2  x2<t) 

=  -x2(t)  -  x^(t) 

_2 _ _ 

s-1 

es 

ES  +  1 


(3.8) 


(3.9) 

(3.10) 


In  this  case  for  which  G  and  AG  are  scalars,  the  condition  (3.7)  is 
equivalent  to 

| 1  +  G~1( j«) |  >  [G( ju,£) i  (3.11) 

or 

|l  +  G(ja»|  >  |G( jw)AG( ja>,e)  ,  (3.12) 

for  all  w>0. 

The  condition  (3.1 2)  has  an  interesting  graphical  interpretation. 
Specifically  the  Nyquist  locus  of  G(j'-j)  must  avoid  the  critical  point  -1 
by  at  least  the  distance  j  G(  joj)AG(  ju) )  J  (Figure  5).  It  is  easily  verified 
that  for  £  «+l  and  <a  »  v2,  we  have  jl  +  G(jco)!  *  |  G(  jcu) AG(  jco)  j  so  that 


A 


-2' 
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Figure  5.  Illustrating  the  Condition  (3.12). 


Example  2 


X1  =  -X1  +  X2 


-x2  =  -  xL  -  x2 


G(s)  = 


(3.14) 


AG(s,e)  = 


(3.15) 


As  in  the  previous  example,  we  will  check  (3.12).  We  have 


/2+u2  >  2  >  1  > 


>  /“ 


so  that  the  system  is  stable  for  all  £  >  01 


(3.16) 


Summarv  and  Conclusions 


In  this  n-aner  we  nave  considered  the  robustness  of  the  stability 
property  of  a  linear  feedback  system  to  variations  in  the  system  model. 

Che  approach  was  to  generalize  a  fundamental  result  in  numerical  linear 
algebra  concerning  robust  inversion  of  matrices  to  linear  operators  of 
the  type  arising  in  input-output  stability  theory.  In  the  time-invariant 
case,  this  approach  specializes  to  give  sufficient  conditions  for  stability 
under  additive  and  multiplicative  perturbations  that  are  easily  verified 
oy  computing  the  singular  values  of  certain  transfer  function  matrices. 

We  then  applied  the  robustness  condition  to  the  analysis  of  singularly 
perturbed  systems.  We  were  able  to  give  an  explicit,  readily  computable 
bound  on  the  magnitude  of  the  perturbation  parameter  £  that  can  be  tolerated 
and  still  have  a  stability  analysis  of  the  reduced  system  valid  for  the  full 
system. 

The  results  of  this  paper  are  felt  to  be  of  interest  for  two  reasons. 
First,  as  Doyle  has  previously  pointed  out  in  the  time  invariant  case  [7], 
the  characterization  and  design  of  robustly  stable  MIMO  feedback  systems 
is  a  fundamental  problem  in  control  theory  that  has  yet  to  be  completely 
resolved.  Second,  as  has  been  previously  emphasized  by  Zames  [18]  and  Safonov 
[5] ,  a  fundamental  problem  in  large  scale  system  theory  is  to  give  conditions 
for  the  success  of  designs  based  on  multiple,  aggregate  models  of  a  single 
large  system  -  this  is  essentially  a  robustness  problem. 
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